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Abstract 

The type of mesh termination used in the context of finite element 
formulations plays a major role on the efficiency and accuracy of the 
field solution. In this work, we evaluate the performance of an ab- 
sorbing boundary condition (ABC) and an artificial absorber (a new 
concept) for terminating the finite element mesh. This analysis is done 
in connection with the problem of scattering by a finite slot array in 
a thick ground plane. The two approximate mesh truncation schemes 
are compared with the exact finite element-boundary integral (FEM- 
BI) method in terms of accuracy and efficiency. It is demonstrated 
that both approximate truncation schemes yield reasonably accurate 
results even when the mesh is extended only 0.3 wavelengths away 
from the array aperture. However, the artificial absorber termination 
method leads to a substantially more efficient solution. Moreover, it 
is shown that the FEM-BI method remains quite competitive with 
the FEM-artificial absorber method when the FFT is used for com- 
puting the matrix-vector products in the iterative solution algorithm. 
These conclusions are indeed surprising and of major importance in 
electromagnetic simulations based on the finite element method. 
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1 Introduction 


The finite element method (FEM) is quite attractive for scattering and ra- 
diation analysis of complex, large scale structures because of its geometrical 
adaptability and low memory requirements. As is the case with any partial 
differential equation method, a difficulty with the FEM stems from the mesh 
termination condition which must be imposed at some distance from the scat- 
terer or radiator. Typically, an absorbing boundary condition is employed 
to terminate the mesh, which is simply an approximate relation between 
the tangential magnetic (electric) field and the tangential derivatives of the 
electric (magnetic) field at the mesh termination boundary. A second order 
ABC retains up to second order tangential field derivatives, the highest to 
be accommodated in an FEM implementation with linear elements. The 
ABC mesh truncation scheme is quite attractive because it leads to sparse 
matrices but, unfortunately, it enlarges the computational domain since the 
ABC must be enforced at some distance away from the surface of the struc- 
ture. The accuracy of the FEM solution improves as the ABC is enforced 
further away from the structure’s surface but, to date, no hard recommenda- 
tions are available as to the optimum/acceptable mesh truncation distance. 
This is because a strong correlation exists among the mesh truncation dis- 
tance, the solution accuracy and the structure’s geometry and composition. 
Among the different 3D FEM implementations, Mayergoyz and D’Angelo 
[1991] employed the Engquist and Majda [1977] ABCs and truncated the 
mesh about one wavelength away from the scatterer. Recently, Chatterjee, 
etc. [1993] showed that the ABC given by Webb and Kanellopoulos [1989] 
can be enforced 0.3 wavelengths away from the scatterer with acceptable so- 
lution accuracy. This ABC also leads to symmetric FE matrices and will be 
employed in this study. 

Alternatively, the boundary integral method can be used for truncating 
the mesh leading to the finite element - boundary integral method (FEM- 
BI). The latter is a numerically exact method but has the disadvantage of 
rendering a partly full, partly sparse system whose solution may require 
greater CPU time and storage. However, if the BI is enforced on special 
mesh termination surfaces (i.e., cylindrical or flat), the fast Fourier transform 
(FFT) can be used in conjunction with an iterative solution of the FEM- 
BI system to maintain a low O(N) storage and minimize the CPU time 
requirements. This has been demonstrated by Jin and Volakis [1991a, b] for 
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the analysis of cavity backed apertures and slots. 

Recently, a third mesh termination scheme was proposed by Jin etc. 
[1992] for two-dimensional scattering. This termination, too, is not exact 
and relies on the use of a metal backed material absorber to terminate the 
mesh. Because this absorber can not, generally, be manufactured, it is often 
referred to as an artificial absorber. In this paper, we propose a metal backed 
artificial absorber (AA) for terminating 3D finite element meshes, and exam- 
ine its accuracy. A comparative study on the computational efficiency of the 
FEM-ABC, FEM-BI and FEM-AA methods is also performed for scattering 
by 3D finite slot arrays. The results of this study are indeed surprising and 
of major importance to the electromagnetics community making use of the 
finite element methods for scattering and radiation. 


2 Geometry Description 

The geometry which was used to perform the comparative study on the effi- 
ciency of the FEM-BI, FEM-ABC and FEM-AA analysis methods is shown 
in Figure 1. It is a finite slot array in an otherwise thick, metallic ground 
plane. The slots are of thickness d and, for this study, they are assumed 
to have a square aperture of width W. There are, of course, no restrictions 
placed on the dimensions and geometry of the slot aperture by any of the 
above methods. Such restrictions result only from the choice of the vol- 
ume elements used in the discretization of the slot volume. For rectangular 
shaped volumes, the simple rectangular elements are sufficient to retain geo- 
metrical fidelity and will be used in connection with all three finite element 
formulations employed in this study. Specifically, the edge-based finite el- 
ement formulation is employed. This formulation is described in [Jin and 
Volakis, 1991b] and explicit expressions for the matrix elements are given 
in [Volakis, etc., 1994]. We also remark that the associated computer codes 
were written such that, the boundary conditions on all metallic surfaces and 
sections separating the slots are enforced apriori before construction of the 
final FE system. Below, we describe other details specific to each of the three 
formulations. 
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3 Description of the Analysis Methods 

Figure 2 displays cross sectional cuts of the 3D slot array and illustrates the 
location of the mesh termination boundary used in connection with the appli- 
cation of the FEM-BI, FEM-ABC and FEM-AA implementations. Through- 
out this study, it will be assumed that the array is illuminated by a plane 
wave incoming from the upper half space. Given such an excitation, of inter- 
est is the evaluation of the reflected and transmitted fields. For the FEM-BI 
implementation, the boundary integral termination is enforced on the top 
and bottom apertures of each slot comprising the finite array. The specifics 
of the implementation are presented in [Jin and Volakis, 1991a,b], and the 
resulting system is solved via the biconjugate gradient method (BiCG). To 
retain an 0(N) memory requirement, the fast Fourier transform (FFT) was 
employed in the BiCG algorithm to perform the matrix- vector products asso- 
ciated with the BI subsystem. The specifics of this implementation are out- 
lined in [Jin and Volakis, 1992] and we note that use of the FFT avoids the 
generation of the full BI subsystem altogether, except for the non-redundant 
matrix elements. It also leads to substantial speed-up in the execution of the 
matrix- vector products since the operations are reduced from O(Njf) down 
to 0(Nb\nNb), where Nb is the number of unknowns on all slot apertures. 

For the implementation of the FEM-ABC formulation, the mesh at the 
bottom slot apertures was extended 0.3 wavelengths beyond the slot array 
surface as illustrated in Figure 2(b). The mesh was then truncated by im- 
posing the second order ABC 

ixH = -y.E, + ^Vx[i(VxE).l + ^V,(V.*,) (1) 

whose accuracy was already examined in [Chatterjee, etc., 1993]. Gener- 
ally, three brick elements were used to fill the mesh between the slot array 
aperture and the ABC surface. As noted on Figure 2(b), the upper mesh 
surface was still terminated by the boundary integral method. This mixed 
truncation scheme was selected to minimize the differences among the three 
implementations, thus, allowing a somewhat more controlled comparison. 
The implementation of the ABC truncation scheme on the upper half space 
would require special considerations to remove the large specular field con- 
tributions from the ground plane. Since this study is concerned with the 
relative merits of the three different termination schemes, it was not deemed 
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necessary to apply the ABC truncation scheme on the upper half space as 
well. For the solution of the resulting system, the BiCG algorithm was again 
employed. As discussed earlier, the matrix-vector product for the BI sub- 
system pertaining to the upper aperture elements was executed using the 

FFT. 

Figure 2(c) illustrates the mesh termination at the lower half space us- 
ing a thin 0.15A o artificial absorber. The absorber dielectric constants were 
found by minimizing the sum of the TE and TM reflection coefficient mag- 
nitudes of the planar metal-backed layer over all real angles of incidence. 
The Simplex method was employed for this minimization. However, given 
that the reflection coefficients must be minimized over a path of complex 
angles to be effective in absorbing diffracted waves, the minimization only 
over real angles does not necessarily lead to an optimum absorber. Thus, it 
was necessary to consider many different metal-backed absorbers (including 
multiple layer absorbers). The effectiveness of each absorber candidate was 
tested by using it to terminate the mesh at the bottom aperture of a single 
slot. Also, of considerable importance was the selection of a metal-backed 
absorber whose refractive index was sufficiently small to avoid a requirement 
for higher density gridding in the material. Note that the top surface of 
the absorber layer was placed only 0.15A o away from the lower slot array 
aperture. Typically, in our implementation, two brick elements were placed 
between the slot aperture and the layer’s surface. Also, two brick elements 
were used to model the lossy absorber. Note that the impedance of the ab- 
sorber material is equal to that of free space since the relative permittivity 
and permeability are equal. Thus, the absorber simply acts as an attenuator 
of all waves impinging upon its surface. 


4 Comparison of the Mesh Termination Schemes 

Before performing a comparison on the efficiency of the aforementioned mesh 
termination schemes, it is first instructive to establish the accuracy of the 
FEM-ABC and FEM-AA solutions. For this purpose, numerous backscatter 
and bistatic patterns were computed for different slot and array configura- 
tions. In all cases, the truncation schemes illustrated in Figures 2(b) and 2(c) 
yielded results which were in good agreement with the rigorous FEM-B1 so- 
lution. Some example bistatic patterns are given in Figures 3 and 4. Figure 3 
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displays the <Jge and cr^ bistatic radar cross-section (RCS) of a 2x2 array. 
The array is illuminated by a plane wave incident along the direction ( <f>' nc = 
0°, 9 ,nc = 80°) for the age computation and along (<^ mc = 0°, 0' nc = 30°) for 
the ( 744 , pattern. The displayed pattern is taken in the upper xz plane and the 
angle a is measured from the x-axis. Figure 4 displays the a 6e bistatic pattern 
of a 15x2 slot array illuminated by a plane wave at ( <j)' nc = 30°, 6 ,nc = 40°). 
The pattern is taken in the upper <f> = 30° plane and a is measured from 
6 = 90°. For all calculations shown in Figures 3 and 4, the slot size was 
1A o x1A„x0.225A o as shown in Figure 1, and the period was 1.15A 0 in the 
x-direction and 1.225A 0 in the y-direction. 

We observe from Figures 3 and 4 that all three solutions (FEM-BI, FEM- 
ABC and FEM-AA) are in good agreement with each other. As expected, the 
FEM-AA formulation required less unknowns than the FEM-ABC formula- 
tion because all degrees of freedom on the metal backing the absorber are 
set to zero apriori. The number of unknowns used in connection with each 
formulation (and the CPU time required to compute the bistatic pattern) 
are given to the right of the legend on each figure. Of course, the FEM- 
BI formulation required the least unknowns since the mesh is terminated at 
the surface of the bottom aperture. Note that a fourth curve (identified by 
the legend cavity ) is included in Figures 3 and 4, and this corresponds to 
the scattering by the same geometry with the bottom aperture of the slots 
shorted. The inclusion of this curve serves to illustrate the level of power 
which must be absorbed by the artificial absorber or the ABC. 

The number of unknowns used in connection with each of the finite ele- 
ment formulations are given in Figure 5. Not surprising, the unknowns grow 
linearly, but the slope for each of the formulations is different. The FEM-BI 
method has the smallest slope and the FEM-ABC method has the largest. 
The corresponding actual CPU times on an HP9000/750 workstation for cal- 
culating the scattered field (single incidence angle) are given in Figure 6. It 
is seen that the FEM-ABC method demands the most CPU time primarily 
because it is associated with the most unknowns, thus, requiring more iter- 
ations for solution convergence. The substantially better performance of the 
FEM-AA and FEM-BI methods is indeed surprising. In view of the unknown 
comparisons given in Figure 5, the superior CPU time performance of the 
FEM-AA method over the FEM-ABC method stems from improved solu- 
tion convergence and not only from the differences in the degrees of freedom 
between the two methods. 
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The CPU time performance of the FEM-BI is indeed impressive given the 
additional computational burden inherent with this method. Obviously, the 
incorporation of the FFT for computing the matrix-vector products in the 
BiCG algorithm has played a major role in reducing the CPU time. Most 
importantly, the improved performance of the FEM-BI method remains unal- 
tered as the array size increases. For this application, both the FEM-BI and 
FEM-AA methods have comparable CPU performances. But the impressive 
performance of the FEM-BI method can only be achieved for planar mesh 
terminations and provided uniform gridding is employed. However, we can 
contend that the artificial absorber termination can be generalized to termi- 
nate meshes on non-planar surfaces making it quite attractive for a variety 
of applications. 


5 Conclusions 

The presented efficiency comparison of the three mesh termination schemes 
has demonstrated that, although both ABC and Artificial Absorber termi- 
nations had good and comparable accuracies, the Artificial Absorber termi- 
nation was the one challenging the BI termination in terms of CPU time 
efficiency. Therefore, the conclusion is that the Artificial Absorber termi- 
nation is a viable option when the boundary integral termination can not 
be used. Given its rigor, the boundary integral termination is, of course, 
the best choice for the analysis of planar array. This is true, provided the 
FFT is employed for computing the matrix-vector products appearing in the 
iterative solution algorithm. 
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_ -gurc x \ A finite rectangular slot array in a thick metallic ground plane 
fa) three dimensional perspective view 
(b) side view of an isolated array cell 
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Figure 2: Illustration of the mesh termination schemes used in connection 
with each of the three finite element formulations 
fa) Boundary Integral Termination (FEM-BI) 

(b) ABC Termination (FEM-ABC) 

(c) Artificial Absorber Termination (FEM-AA) 
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Figure 3. Bistatic RCS of a 2x2 slot array whose element geometry and 
separation distance are given in Figure 1. The excitation is 
a plane wave and the pattern is taken in the xz plane 

(a) gqq RCS with (<^> ,nc = 0°, 6 inc = 80°) 

(b) o H RCS with ( <f) inc = 0°, 0 inc = 30°) 
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Figure 4: Bistatic a 0 g RCS of the array shown in Figure 1 in the <f> = 30° 
cut for a plane wave incidence at ( (f> inc = 30°, Q inc — 40°) 
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gure 5. Comparison of the number of unknowns for the three 
termination schemes vs slot array size (see Figure 1) 
tor this comparison, the array has N x x2 elements 
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Figure 6: Comparison of the CPU time required by each of the finite 
element formulations to compute the scattering for one 
incidence angle vs array size (see Figure l).The array is 
shown in Figure 1 and has N x x2 elements 



